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A method is presented for the computation of the far field radiation 
patterns of paraboloid reflector antennas by using a modified stationary 
phase approximation to eliminate one integration. This method is applicable 
to open cassegrain, offset paraboloid and horn reflector antennas. For 
symmetrical paraboloid antennas the modified approximation reduces to 
the exact expression obtained by direct integration. 

The errors introduced by the stationary phase and modified stationary 
phase approximations are investigated. Specifically the far field of an open 
cassegrain with a 128 ivavelength aperture diameter is computed by the 
approximate method up to 20 degrees off-axis. The difference between these 
radiation patterns and those computed by double integration, is less than a 
few tenths of a dB up to 1.0 degree, and less than a few hundredths of a dB 
at larger angles off-axis. 

In order to estimate the computational advantage of this approximation, 
the number of points required for integration of an oscillatory function by 
Simpson's rule is also examined and it is determined that at least 6 points 
per cycle are necessary to obtain 4 decimal accuracy. For fewer points the 
error is appreciable. 

I. INTRODUCTION 

The computation of the far field radiation patterns of large reflector 
antennas is of importance in predicting the performance of satellite 
ground stations. For example, the antenna sidelobes contribute to the 
system noise temperature and may cause interference with other com- 
munication systems. The open cassegrain antenna 1 is a particularly 
suitable configuration for obtaining low sidelobe levels, since blocking 
by the subreflector and its supports are eliminated. A further advantage, 
resulting from this feature, is that the radiation pattern can be accu- 
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rately predicted and it has been shown to be in good agreement with 
experimental results. 1 

The previous computations of the radiation patterns of open cas- 
segrain antennas have been performed by precise computation of the 
appropriate diffraction integrals, generally requiring a double numerical 
integration. For large angles off-axis, such computations require con- 
siderable computation time. It is, however, for large angles that the 
integrals which are used for the computation of the far field radiation 
patterns are of a form which is suitable for approximation by the method 
of stationary phase. This method was initially applied to eliminate the 
azimuthal <t> integration, but it was subsequently recognized that certain 
terms in the approximation are related to the asymptotic expansions of 
Bessel functions. The stationary phase approximation could therefore 
be modified, with the observed result that the far field radiation pattern 
can be computed with good accuracy also in the immediate vicinity of the 
main beam. 

In the following sections we derive the stationary phase approximation 
and present a geometrical interpretation of the location of the station- 
ary points. For the far field on axis, the tf> integration is performed in 
closed form and it is shown that the antenna gain is the same for both 
perpendicular polarizations. Numerical computations are performed to 
estimate the error introduced by the stationary phase and modified 
methods. The extended range of applicability of the latter method is 
evident from the computations. The number of points per cycle needed 
to obtain an accurate value for an integral of an oscillatory function is 
also examined, and it is shown that for the functions considered at 
least 6 points/cycle are needed. 

The far field radiation patterns of an open cassegrain with a 128 
wavelength aperture diameter are computed with this method up to 20 
degrees off-axis. In the vicinity of 20 degrees, the relative sidelobe levels 
are less than — 65 dB or about 15 dB below isotropic. 

l.i The Far Field 

The far electric field E, of a paraboloid reflector antenna in an angular 
region about the axis, can be, based on the projected aperture field 
method, related to the reflected field at the aperture, E r , by the follow- 
ing expression: 2 

E ' = ieXP x£ j7bfin) IL Er( *" ' Vv) exp (jk9 "' ljJ ds (1) 
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where 



X = wavelength, 

k = 2tt/X (propagation constant), 
R a = the distance to the far field observation point, 
q p is a vector in the aperture plane, 
l Ra is a unit vector which specifies the direction of the observation 
point, 
A is the aperture area. 

Specifically the direction of the observation point \ Ra expressed in 
terms of the unit vectors of the aperture (x p , y p , z p ) coordinate 
system is: 

1«. = 1* P sin 6 a cos <t> a + I,, sin O sin <f> a + 1., cos d a (2) 

where 6 a and <j> a are the far field observation angles, and 

p p = l,^Cp + KVv ■ W 

For an open cassegrain the incident fields at the main reflector can 
be more readily computed in a spherical coordinate system with the 
axis aligned with the horn subreflector axis as shown in Fig. 1. There- 
fore, the integrations in equation (1) are also performed in this co- 
ordinate system. The relations between aperture coordinates and the 
fields in the two coordinate systems were derived previously 1 and are 

x„ = r[cos 6 sin 6 cos <t> + sin 6 cos 6], (4) 

y p = r sin 6 sin <f>, w 

*L E r = 1 { [sin d Q sin d - cos 4>(1 + cos 6 cos o )]E e 
r 

+ sin 0(cos 6 + cos )E+} 

- I„ p {sin0(cos 6 + cos B)E e 

- [sin 0sin 0„ - cos^(l + cos d cos 9 )]E^}, (6) 

where E„ and E+ are the 6 and <t> components of the incident electric 
field, / is the focal length of the paraboloid and r is the equation of the 
paraboloid surface in the d, <f> coordinate system 

. m 2/ ; (7) 

1 + cos do cos 6 — sin 6 sin 6 cos </> 
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Fig. 1 — Open cassegrain antenna. 

O is the offset angle. The surface element 

ds = r 2 sin 6 dd d<f>. (8) 

Heretofore the above integral has been evaluated by double inte- 
gration using Simpson's rule. Although rather accurate results can be 
obtained in this fashion, computation time for a given angle, 6 a in- 
creases roughly proportional to (sin 0„) 2 (see Appendix C). As a result, 
except for the mainlobe and first few sidelobes of the far field, this ap- 



STATIONARY PHASE METHOD 435 

proach becomes extremely time consuming. To reduce the computation 
time for large values of the off-axis angle 0„ , the stationary phase 
approximation to the <f> integration is investigated. 

1.2 The Method of Stationary Phase 
Consider an integral of the form 

/(A: sin 9„) = I g(4>) exp [j(k sin d a )^)) d<t> (9) 

where (k sin 0„) is large, tK#) is a real function and g(<t>) is a slowly 
varying function. The method of stationary phase 3 approximates the 
above integral to 0(l/k sin a ) by considering only contributions in the 
vicinity of the stationary points <f> t where ^'(0.) = 0. 
Under these conditions 

I(k sin O « E g(0,) ( fcsin 2 ^, (0o y exp {jfe sin *.*(*,)}. (10) 

This method has been apphed to the </> integration in the expression for 
the far field (1). 

From equations (1), (2) and (3) ^(0) can be written as 

f(ft>) = [la. cos a + l„ p sin <£„]■ p p . (11) 

For specified observation angles 0„ and <t> a , equation (11) can be con- 
sidered as the projection of the vector p p in the direction of the unit 
vector l Pa = l Ip cos <p a + l Vp sm<f> a . For the problem under consideration 
p p is a function of 6 and <f>. It has been shown previously 1 that for con- 
stant 6, the vector q p describes a circle as <£ varies from to 2ir. The 
equation of the circle is 

2/ sin O V ■ = / 2/ sin V Q2) 

Xv cos O + cos 0/ ^ y " Vcos + cos 0o/ ' 

A family of such circles for an offset angle O of 47.5° is shown in Fig. 2. 

Therefore the condition t//(tf>) = corresponds to determining the 
extreme values of the projections of the vector q v in the direction of the 
unit vector l Po . It is evident from Fig. 2 that as p„ describes a constant 
circle two extreme values for the projections exist, namely at those two 
points on the circle such that tangents to the circle passing through the 
points intersect normally a line in the direction of the unit vector l Pa . 
Furthermore the difference between the two extreme projections is the 
circle diameter. 

The expressions for the stationary points and the other values which 
enter in the evaluation of the <t> integration are derived in Appendix A. 
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Fig. 2 — Projection circles of the paraboloid reflector for e = 47.5°. 

1.3 Approximate Values 

In order to approximate the <$> integration in equation (1), namely 



E = / r 2 E r exp {jk sin 6 a {x P cos <f> a + y P sin <£ a )} d<}> 

Jo 



(13) 



by the stationary phase method, the reflected field E r has to be de- 
termined at the stationary points. This field has an explicit <f> dependence 
for a TE tl mode or combined TE lx — TM lx excitation. For these modes 
it has been shown 1 that for x and y polarization the field components 
E g and E t are 

, polarization \ E > = **> C0S * CM) 

[E* = -E(0) sin 



y polarization J e 

UEL = E(0) cos <f> 



(15) 
(16) 

(17) 



where E(t/2) and E(0) denote the d dependence of the fields when 
the feed horn is excited for y polarization in the planes <f> = ir/2 and 
<j> = 0, respectively. 
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As shown in Appendix A, it is sufficient to evaluate the fields at the 
stationary points for one polarization only. For the other polarization 
the fields are then readily obtained. Therefore only y polarization will 
be considered. Furthermore it will be assumed that the radiated field 
from the subreflector has been computed at a constant radius from the 
focal point of the hyperboloid subreflector. 

By assuming a 1/r dependence for subreflector fields, the relation 
between the fields is 



E r (0, 0) = 



2/ E r (0, <t>) 



1 + cos 



(18) 



where E c is the field at the distance 2//(l + cos O ) and r is the equation 
of the paraboloid (7). The spherical phase dependence of the field in 
equation (18) is suppressed. 

With the approximation (18) and the stationary phase approximation 
to the <f> integration, the far field using equation (1) is obtained from 
the integral (see Appendix A) 



(B / .)-i exp( -f R - ) / , "(E,)8in^g 



(19) 



where the subscript y designates that the far field is for y polarization, 
and 8 m is the illumination angle. The y and x (cross polarization) com- 
ponents of (E„) namely (E v ) v and (E v ) x are given by 



W, = 



-T(2/)y 



c(a 2 — b 2 cos 2 <£„)(! + cos O ) 



a + 4/J 



• [c(e ia + e- ia ){a sin 2 0„ E e (*/2) + c cos 2 <f> a E e (0) J 
- &sin 2 0„ cos0 a (e' a - e~ ia ){cE c (^/2) - aE c (0)}] 
and the cross polarized component (E v ) x is 



(20) 



(E v ) x = 



ir(2/)V»sin6, 



c(a — b cos </>„)(! -(- cos a ) 



""l« + ^yj 



where 



[b(e ia ~ e~ ia ){a sin 2 O E e (0) + c cos 2 O E c (tt/2)\ 

+ c cos<t> a (e ia + e- ia ){cE c (0) - aE e (W2)}] 

a = 1 + cos 6 cos do, 



(21) 
(22a) 
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b = sin 0sin 9 , (22b) 

c = cos 6 + cos 0o- (22c) 

2kj sin 0„ sin ,. /onN 

« = — n r ~~ T / 4 - (w 

cos + cos O 

P = 24/ sin e. cos * .[— JfeA_ ™J- - 1 • (24) 

' Lcos + cos O cos O + cos m J 

The second term in equation (24) reduces the phase by a constant. 
The expansions (20) and (21) are valid for (k sin a ) very large. It 
should be noted that (20) and (21) display singularities at a = — ir/4, 
that is, at O = 0, or = 0. However, upon examination it can be seen 
that (20) and (21) contain the first terms of the asymptotic expansions 
for Bessel functions, that is, 

Ux) ~ w cos ( x - f - ;) (25) 

~ i w "(^) \ ( - ir exp X s " *). + exp M' _ 5 

By identifying and replacing the asymptotic terms by the actual 
Bessel functions, the singularities are removed and it might be ex- 
pected that the approximation for small and a would improve and 
furthermore for large values of and 9 a the approximations would be 
equivalent. There is however, no unique method to introduce such a 
replacement. The method chosen was dictated by the requirement that 
for the symmetric case O = 0, the expressions reduce to the exact 
expressions previously determined. 4 This necessitates associating a 
Bessel function of order n, J n (x) with terms cos n^> a or sin ruf>., . 

On this basis the approximations to the <£ integration are: 

(F) -T(2flV' 

{ v,v c(a 2 - b 2 cos 2 0„)(1 + cos O ) 

•Wo{x)[aE t (T/2)+cEM\ +/ 2 (z)cos24 a {a£ c (7r/2) ~cE c (0)}) 
- j/2b{cE e (ir/2) -aE c (0)}{Jx(x) cos <f> a + J 3 (x) cos 30 a }] (26) 

and the cross polarized component 



(EX = 



c(a 2 - b cos O ) (1 + cos O ) 

• [jb{3aE c (0) + cE e (T/2)}Jy(x) sin 0„ 

+ 2c{aE e (ir/2) - cE c (0)} J&) sin 2tf a 

+ jb{aE c (0) - cE e (ir/2) } J 3 (x) sin 30 Q ] (27) 
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where 

_ 4tt/ / sin fl„ sin 6 \ ^ 

X Vcos + cos 6 J 

It is noted the (E y ) y is symmetrical with respect to <f> a since by inter- 
changing 4> a by — O the expression remains the same. This would be 
expected since the plane <f> = is a plane of antenna symmetry. It is 
also evident from equation (27) that the cross polarized component is 
zero in the plane of symmetry and is antisymmetrical with respect to O . 

The above expressions reduce to those obtained by the method of 
stationary phase for large values of x. As shown subsequently by nu- 
merical integration the latter approximations extend considerably the 
range of and a beyond which the stationary phase approximations are 
applicable. 

For x polarization as outlined in Appendix A the expressions are 
similar. In particular {E x ) x is of the same form as (E y ) v with E c (o) and 
E e (ir/2) interchanged. The cross polarized component (E x ) y is of the 
same form as —{E y ) x with E c (o) and E c (ir/2) interchanged. 

The expressions (26) and (27) are of course approximate. This is 
evident by considering the special case 0„ = 0, where the values for the 
fields must be the same independent of 4> a . For this special case the <f> 
integration is performed in closed form in Appendix B. 

For the antenna shown in Fig. 1, with O = 55° and m = 34.0°, by 
assuming the radiation fields of the subreflector at a constant distance 
are the same in the E and H planes, it is shown that the differences 
between the exact and the approximate values at O = are 0.049 in 
the .ft-plane and —0.053 in the #-plane, both in comparison to one. 
The subsequent numerical computations indicate that these are the 
largest errors introduced by the approximation. 

1.4 Numerical Results 

In order to determine the validity of the above approximations, 
computations of the far field radiation patterns for the open cassegrain 
antenna have been performed using the subreflector radiation pattern 
E e shown in Fig. 3. 

The integration with respect to 0, as indicated in equation (13), has 
been performed, employing Simpsons rule, as a function of in the E 
and H planes, and for observation angles 6 a = 0, 2.5, 5°, 10° and 20°. 
Estimates for the number of points required for the <j> integration and 
the computation time are presented in Appendix C. 

The normalized amplitudes obtained from the integration are shown 
in the upper portions of Figs. 4-8. The normalization was based on the 
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16 20 24 28 32 36 

ANGLE OFF AXIS IN DEGREES 



Fig. 3 — Amplitude and phase of subreflector radiation pattern. 



stationary phase approximation (20) which shows that in the planes 
O = or 7r/2 the integral (13) is proportional to E e (o) or E c (x/2) re- 
spectively. The normalized values for the integrals E N shown are 



E„ = 



(EX/E C 
[(E u ) v /E c ]o 



(29) 



where the subscript zero indicates the value at = 0. 

Immediately beneath E N in Figs. 5-8 is shown a plot of the absolute 
value of the difference between the normalized values obtained by inte- 
gration and the stationary phase approximation given by equation 
(20). The third plot in each figure show the corresponding difference 
using the modified stationary phase approximation (26). As predicted 
by the method of stationary phase and as shown in Figs. 5-8, the 
approximations improve as k sin d a increases. However, where as 
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for small values of 6 the stationary phase approximation (20) introduces 
significant error, the modified approximation (26) becomes increasingly 
accurate. It should also be noted that with both approximations rela- 
tive maximum differences occur near zeros, therefore resulting in less 
significance in the second integration. 

Figures 9 and 10 show the amplitudes of the far field radiation as 
computed by single integration and the modified stationary phase 
approximation (26). Shown for y polarization are the far field in the 
plane of antenna symmetry <p a = and jr, and the fields in the plane of 
asymmetry <jy a = ir/2 up to 20° off axis. Figure 11 shows the difference 
in the plane <£ a = between the far field pattern computed by the ap- 
proximate method and the same pattern computed using double integra- 
tion. Excluding the vicinity of relative minima, errors were less than 0.2 
dB up to off-axis angles of 1°, and on the order of a few hundredths of a 
dB for larger angles. It should be noted that on axis the difference is zero, 
since the exact expression for the <f> integration as given in Appendix B 
is incorporated in the single integration program. 

II. CONCLUSIONS 

A method has been developed for the numerical computation of the 
far field radiation patterns of open cassegrain antennas and related 
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Fig. 4 — Error introduced by modified approximation [equation (26)] on axis. 
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Fig. 5 — Error introduced by stationary phase approximation [equation (20)] 
and modified approximation [equation (26) ] . 

antenna configurations using a single numerical integration, which re- 
duces considerably the computation time. The method is based on the 
stationary phase approximation but modified such that for symmetrical 
paraboloid antennas the approximation reduces to the exact expression 
which is obtained by direct integration. The errors introduced by 
stationary phase and modified approximations are examined. It is shown 
by numerical computations that the error in the far field radiation 
pattern introduced by the modified stationary phase approximation at 
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angles beyond the main beam is on the order of a few hundreths of a dB. 

The number of points needed for numerical integration of oscillatory 
functions by Simpson's rule is also examined by using specific oscillatory 
functions. It was determined that at least 6 points per cycle are necessary 
to obtain four decimal-point accuracy. 

The radiation patterns of an open cassegrain antenna with a 128 
wavelength aperture diameter are computed up to 20° off axis. In this 
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Fig. 6 — Error introduced by stationary phase approximation [equation (20)] 
and modified approximation [equation (26)] 
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Fig. 7 — Error introduced by stationary phase approximation [equation (20)] 
and modified approximation [equation (26)]. 



region the sidelobe levels are —65 dB or —15 dB below the isotropic 
levels. 

Although the stationary phase method has been applied here to the 
computation of the far field based on the projected aperture field method, 
the same approach may be used for the computations based on the 
current distribution method. In particular in the plane of antenna sym- 
metry, the locations of the stationary phase stationary points are the 
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same and therefore the presented approximations are readily extended 
to this plane. 



APPENDIX A 



Derivation oj the Stationary Phase Approximation 
Consider the <£ integration of equation (1) 
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Fig. 8 — Error introduced by stationary phase approximation [equation (20)] 
and modified approximation [equation (26)]. 



446 



THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1970 



«2* 

E(k, 0) = / E r (0, 4>) exp [jk sin 6 a (x„ cos <t> a + y p sin </> a )> 2 d<t> (30) 

with x v ,y v , E r and r given by equations (4) through (7) . The stationary 
points are determined by 



— cos* + — sin* a = ^ = 0. 



(31) 



This leads to the relation 

[(1 + cos cos d ) cos <f> — sin sin O ] sin tf> 

— sin <£[cos 6 + cos O ] cos a = 0. (32) 
Equation (32) is a quadratic equation in cos <f>. Solving this equation 
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Fig. 9 — Far field radiation pattern in the plane of symmetry. 
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-16 




19.5 21.0 



Fig. 10 — Far field radiation pattern in the plane of asymmetry. 



yields the following values for the stationary points 

_ sin sin n ± cos <£„(! + cos cos O ) 
cos li8 -= l + cQs e cog ^ ± cog ^ gin e gin do 



and from (32) 



sin 0, i2 = =1= 



sin <ft n (cos + cos O ) 

1 + cos 8 cos O ± cos a sin 8 sin O 



Evaluating the phase factor \f/ gives 

, . . n , (si n 8 + cos a sin O ) 
* (</,l) = 2/ COS 0+ COS 0o "' 

(sin — cos <ft a sin O ) 
^ 2) = ~ 2/ COS + COS O 



(33) 
(34) 

(35) 
(36) 
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Fig. 11 — Difference in far field pattern as computed by double integration and 
approximate method (in the plane #« = 0). 



Note that 



™-*>-telt2ni 



(37) 
.cos d + cos O ' 

which is the diameter of the circle given by equation (12). That is the 
stationary points are antipodes on the projection of the plane of the 
intersection of the paraboloid surface with the cone = constant. 
Evaluation of the second derivative leads to 

[1 + cos cos O ± cos <j> n sin sin d f 






= =F2/sin 



[cos 4- cos d ]' 



Evaluation of r gives 



ri.2 = 2/ 



[1 + cos cos dp dr sin sin O cos <fr n ] 
[cos 6 + cos 6 ] 2 



(38) 



(39) 



It remains to evaluate E r at the stationary points. From equation (18) 

E, 2/ 



E r = 



(40) 



r (1 + cos 6 ) 
From equations (6) and (31) 

(jy, , _ _S*4m feg * + «» «J [!„(£«, cos *. - E rt sin *.) 
v r/1,a sin0 a 1 + cos d 

+ 1 W (E ## sin </,„ + # e , cos 4>,)] ■ (41) 
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Substituting the values for x and y polarization in equations (14) 
through (17) it is evident that the expressions are similar. Therefore 
the x component for x polarization can be obtained from the y com- 
ponent for y polarization by interchanging E e {o) with E(ir/2). Similarly 
the y component for x polarization can be obtained from the x com- 
ponent for y polarization also by interchanging E c (o) with E(ir/2) and 
changing the sign in front of the resulting expression. 

Substituting the appropriate expressions for y polarization into the 
stationary phase approximation leads to equations (20) and (21). 

APPENDIX B 

The <j> Integration on Axis (0 O = 0) 

From equations (6), and (13) for d a = 0, the integral for the y com- 
ponent of y polarization with respect to <t> can be written 

w« (i + C o S e ) 

f 2 * i \E r (r/2)c - aEM] sin 2 <i> , EM \ ^ ( 42 ) 

J \ [a - b cos </>] 2 a - b cos <*>/ 

where a, b, and c are denned by equations (22a), (22b), and (22c). 
Integrating the first term by parts, reduces the evaluation of equation 
(42) to a tabulated integral, 5 that is, 

— T-. i = r^w < 43) 

a — b cos <f> (a — b )' 

(F \ - 2(2tf {EM2) + EM] (44) 

(E r ,) y - -2(21) (1 + cQs ^ )2(1 + cog Q) (44) 

The integration for x polarization (E TZ ) X gives the same result. As a 
consequence the on-axis gain for an open cassegrain antenna is the same 
for x and y polarization if the excitation is the same. Based on the 
approximation (26) 

(2/) 2 [ff,W2)(l + cos g cos go) + E r (0)(cos g + cos fl„)] .„. 
{ v) » (1 + cos O )[(1 + cos 6 cos O ) 2 - (sin 0sin d Q cos <£„) 2 ] 

To estimate the relative error it is assumed that E c (o) = E e (ir/2). 
This gives in the plane <f> a = or ir 

[ (I + cos gXl + cos O ) 1 2 



hence 
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Fig. 12 — Numerical integration of oscillatory function. 



Similarly in the plane <j> a = ir/2 

(AF \ _ i _ [~ (1 + cos g)(l + cos d ) 
K v)v 2(1 + cos 6 cos e Q ) 

Equations (48) and (49) can be combined yielding 

1 



]•• 



(All.). = 1 - 



[1 ± (tan 9/2 tan O /2)T 



(46) 



(47) 



where the minus sign corresponds to <f> = or ir and the plus sign is for 
0. = t/2- 

appendix c 

Computation Time Estimates 

The attempt to find a suitable approximation to be used for the 
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evaluation of far field radiation patterns was motivated, in fact necessi- 
tated, by the excessive computation time required for a double inte- 
gration procedure. 

Both the integral [equation (13)] and the integral [equation (19)] 
have oscillatory integrands, with the <f> integrand, for the maximum 
value of 0, having approximately double the number of oscillations as 
the integrand. The maximum number of full cycles in the integrand 
is equal approximately to 

§sin* a 

where D is the antenna aperture diameter. 

Figures 12 through 15 show the <f> integrand for 6 = 34°, and for 
various values of observation angle O . Below these figures is shown 
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Fig. 13 — Numerical integration of oscillatory function. 
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Fig. 14 — Numerical integration of oscillatory function. 

the result of the <£ integration using Simpson's rule with increasing 
number of points. We have found by numerical integration that for the 
integrands discussed above, a minimum number of 6 points per cycle is 
necessary to provide reasonable accuracy. 

For a discussion of numerical integration of oscillatory integrands 
the reader is referred to Ref. 6. 

An estimate for the processor time required to calculate both com- 
ponents of the far field pattern in both the E and H planes (4 patterns) 
out to an observation angle of d a degrees by the method of double inte- 
gration is given by: 

*A». 3 

where 
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k 2 is the fundamental program loop execution time (1.13 X 10~ 4 

min., on the G. E. 635). 
N is the number of integrand evaluations per cycle required by the 

integration procedure. 
A0„ is the observation angular increment at which results are to be 

calculated. 

An estimate for the same calculation by the approximate method 
derived herein is: 



where 



h —2- N — 



fc, = 0.75 X 10" 3 min. 
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Fig. 15 — Numerical integration of oscillatory function. 
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Table I — Execution Times 



0o 


5° 


10° 


20 


Double Int. 
Appr. Method 


4.7 min. 
0.9 min. 


38 min. 
3.7 min. 


300 min. 
15 min. 



Table I shows comparable execution times for various observation 
angle extremes, assuming Ad a = 0.1°, N = 10. 
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